Electrically detected interferometry of Majorana fermions in a topological insulator 
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We show how a chiral Dirac fermion (a massless electron or hole) can be converted into a pair 
of neutral chiral Majorana fermions (a particle equal to its own antiparticle). These two types 
of fermions exist on the metallic surface of a topological insulator, respectively, at a magnetic 
domain wall and at a magnet-superconductor interface. Interferometry of Majorana fermions is a key 
operation in topological quantum computation, but the detection is problematic since these particles 
have no charge. The Dirac-Majorana converter enables electrical detection of the interferometric 
signal. 
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There is growing experimental evidence [TJ [3] that 
the 5/2 fractional quantum Hall effect (FQHE) is de- 
scribed by the Moore-Read state [4J. This state has re- 
ceived much interest in the context of quantum compu- 
tation [5], because its quasiparticle excitations are Ma- 
jorana bound states. A qubit can be stored nonlocally 
in a pair of widely separated Majorana bound states, so 
that no local source of decoherence can affect it [BJ. The 
state of the qubit can be read out and changed in a fault- 
tolerant way by edge state interferometry [7, 8, 9 . This 
"measurement based topological quantum computation" 
[TU] combines static quasiparticles within the Hall bar to 
store the qubits, with mobile quasiparticles at the edge 
of the Hall bar to perform logical operations by means of 
interferometric measurements . 

The electronic correlations in the Moore-Read state in- 
volve a pairing of spin-polarized fermions, equivalent to a 
superconducting pairing with p x + ip y orbital symmetry 
[TT1 [T2l IT5] , Such an exotic pairing might occur nat- 
urally in the Si'2Ru04 superconductor |14j . or it might 
be produced artificially in p-wave superfluids formed by 
fermionic cold atoms [15]. Recently, Fu and Kane [T6] 
showed how a conventional s-wave superconductor might 
produce Majorana bound states, if brought in proxim- 
ity to a topological insulator. This class of insulators 
has metallic surface states with massless quasiparticles, 
as has been demonstrated in B^Sbi-^ alloys |17j and 
Bi2Se3 single crystals [T5| [I5]. The latter material is par- 
ticularly promising for applications because it remains 
a topological insulator at room temperature. The 5/2 
FQHE, in contrast, persists only at temperatures well 
below IK [DUE]. 

While induced superconductivity in a topological in- 
sulator seems an attractive alternative to the FQHE for 
the purpose of quantum computation, one crucial differ- 
ence creates a major obstacle: Quasiparticle excitations 
carry a charge in the FQHE, but they are charge-neutral 
in a superconductor. All known schemes [3 019] for edge 
state interferometry rely on electrical detection, and this 
seems impossible if the edge states carry no electrical 
current, ft is the purpose of this work to propose a way 
around this obstacle, by showing how a pair of neutral 
Majorana fermions can be converted phase coherently 
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FIG. 1: Three-dimensional topological insulator in proximity 
to ferromagnets with opposite polarization (M-f and M|) and 
to a superconductor (S). The top panel shows a single chiral 
Majorana mode along the edge between superconductor and 
ferromagnet. This mode is charge neutral, so it cannot be 
detected electrically. The Mach-Zehnder interferometer in the 
bottom panel converts a charged current along the domain 
wall into a neutral current along the superconductor (and 
vice versa). This allows for electrical detection of the parity 
of the number of enclosed vortices, as explained in the text. 



and with unit probability into a charged Dirac fermion. 

We first give a qualitative description of the mech- 
anism of electrically detected Majorana interferometry, 
and then present a quantitative theory. Our key idea is 
to combine edge channels of opposite chiralities in a sin- 
gle interferometer, by means of a magnetic domain wall. 
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The appearance of counterpropagating edge channels in 
a single superconducting domain is a special feature of 
a topological insulator in proximity to a ferromagnet, 
where the propagation direction is determined by the way 
time reversal symmetry is broken outside of the conden- 
sate (hence by the polarization of the ferromagnets) - 
rather than being determined by the order parameter of 
the condensate (as in &p x ±ip y superconductor or FQHE 
droplet) . 

Refering to the lower panel of Fig. [I] we see that elec- 
trons or holes (with Dirac fermion operators c a and c a ) 
propagate along the domain wall a until they reach the 
superconductor, where they are split into a pair of Ma- 
jorana fermions 7& and j c of opposite chirality: 



lb + ijc, Ca^lb- ijc- 



(1) 



(Here we have used that 7 = 7^, which is the defining 
property of a Majorana fermion.) 

The Dirac-to-Majorana fermion conversion expressed 
by Eq. relics on the fact that the electron or hole 
mode at the domain wall couples to a pair of Majo- 
rana modes, so that the full information encoded by the 
complex fermion c a is encoded by two real fermions jt, 
and 7c. This is the essential distinction from the pro- 
cess of electron tunneling into a Majorana bound state 
[2171 |2"T1 |2"21 |2"5] . which couples to a single Majorana 
fermion and can therefore not transfer the full informa- 
tion. 

Upon leaving the superconductor the Majorana 
fermions recombine into an electron c* d or hole Cd de- 
pending on the number n v of superconducting vortices 
enclosed by the two arms of the interferometer, 



Cd- 



(2) 



For n v an even integer, no charge is transfered to the su- 
perconductor, while for n v odd a charge ±2e is absorbed 
by the superconducting condensate. The conductance G, 
measured by application of a voltage between a point on 
the domain wall and the superconductor, becomes equal 
(in the zero-temperature, zero- voltage limit) to G = for 
n v = even and G = 2e 2 jh for n v — odd. 

Proceeding now to a theoretical description, we recall 
that the surface of a three-dimensional topological in- 
sulator, in the presence of a magnetization M(r) and 
superconducting order parameter A(r), is described by 
the following Hamiltonian [16 : 
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Here p — (p x ,p y ,0) is the momentum on the surface, 
er = (a Xl Oy, a z ) is the vector of Pauli matrices, vf 
is the Fermi velocity, and Ep the Fermi energy. The 
two magnetizations Mj and Mj in Fig. [T] correspond 
to M = (0,0, M ) and M = (0,0, -M ), respectively. 
Particle-hole symmetry is expressed by the anticommu- 



tation = —z,H of the Hamiltonian with the operator 



-i<jyC 



i<jyC 



(4) 



with C the operator of complex conjuation. 

There is a single chiral Majorana mode with amplitude 
■0 (group velocity v m ) at a boundary between a region 
with a superconducting gap and a region with a magnetic 
gap |16j . At a domain wall between two regions with 
opposite signs of M z there are two chiral Dirac fermion 
modes, an electron mode with amplitude <p e and a hole 
mode with ampitude <f> . The scattering matrix 5i n (e) 
describes scattering at excitation energy e from electron 
and hole modes (along edge a) to two Majorana modes 
(along edges b and c in Fig. [T|, according to 
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Particle-hole symmetry for the scattering matrix is ex- 
pressed by 



1 

1 



(6) 



At small excitation energies |e| <§C |M Z |,|A| the e- 
dependence of S{ n may be neglected. Then Eq. ^ to- 
gether with unitarity (S^ 1 — Sj n ) fully determine the 
scattering matrix, 



5in _ V2 -i 



e la 
e- ia 



(7) 



up to a phase difference a between electron and hole 
(which will drop out of the conductance and need not 
be further specified). 

The scattering matrix S on t for the conversion from Ma- 
jorana modes to electron and hole modes can be obtained 
from S ln by time reversal, 



S out (M) = Sl(-M) 



= J_ fe ia ' \ (1 



V2 V e" lQ / V 1 



(8) 

The phase shift a' may be different from a, because of 
the sign change of M upon time reversal, but it will also 
drop out of the conductance. 

The full scattering matrix S of the Mach-Zehnder in- 
terferometer in Fig. [T] is given by the matrix product 

r< ( S ee S e }i \ ^ f €. ^ h \ n (r\\ 

S = [S he S hh ) ~ bout { J bin > (9) 

where (3b and (3 C are the phase shifts accumulated by the 
Majorana modes along edge b and c, respectively. The 
relative phase 



Pb — C = eSL/hv m + 7r + n v n 



(10) 
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from the transfer matrix 



FIG. 2: Fabry-Perot interferometer, allowing to measure the 
state of a qubit encoded in a pair of vortices. Black lines 
represent electron or hole modes at domain walls, gray lines 
represent Majorana modes at magnet-superconductor inter- 
face. 



consists of three terms: A dynamical phase (proportional 
to the length difference SL = Lb — L c of the two arms of 
the interferometer), a Berry phase of 7r from the rotation 
of the spin-1/2, and an additional phase shift of ir per 
enclosed vortex. 

The differential conductance follows from 



G(V) = ^-\S he (eV)\ 2 =^-sm 2 
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As announced in the introduction, the linear response 
conductance G(0) vanishes if the number of vortices is 
even, while it has the maximal value of 2e 2 /h if the num- 
ber is odd. 

The Mach-Zehndcr interferometer can distinguish be- 
tween an even and an odd number n v of enclosed vortices. 
The next step towards measurement based topological 
quantum computation is to distinguish between an even 
and an odd number rif of enclosed fermions. If n v is 
odd, the parity of n/ is undefined, but if n v is even, the 
parity of n f is a topologically protected quantity that de- 
termines the state of a qubit 5 . To electrically read out 
the state of a qubit encoded in a pair of charge-neutral 
vortices, we combine the Fabry-Perot interferometer of 
the FQHE [5J |9] with our Dirac-Majorana converter. 

The geometry is shown in Fig. [2} Electrons are injected 
in the upper left arm a of the interferometer (biased at 
a voltage V) and the current / is measured in the upper 
right arm e (which is grounded) . The electron at a is split 
into a pair of Majorana fermions ipb and ip c , according to 
the scattering matrix S- ln . A pair of constrictions allows 
tunneling from ij) c to tpd, with amplitude tdc- Finally, 
the Majorana fermions tpd and ipb are recombined into an 
electron or hole at e, according to the scattering matrix 
5 out . The resulting net current I = (e 2 / h)V (\T ee \ 2 - 
\The\ 2 ) (electron current minus hole current) is obtained 
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Notice that the current is proportional to the tunnel am- 
plitude, rather than to the tunnel probability. In the 
low- volt age limit, to which we will restrict ourselves in 
what follows, the phase shift (3b vanishes and tdc is real 
(because of electron-hole symmetry) — so / directly mea- 
sures the tunnel amplitude. 

In general, two types of tunnel processes across a con- 
striction contribute to tdc- A Majorana fermion at the 
edge of the superconductor can tunnel through the super- 
conducting gap to the opposite edge of the constriction 
either directly as a fermion or indirectly via vortex tun- 
neling [24 . Fermion tunneling typically dominates over 
vortex tunneling, although quantum phase slips (and the 
associated vortex tunneling) might become appreciable in 
constrictions with a small capacitance |25j . Only vortex 
tunneling is sensitive to the fermion parity rif, through 
the phase factor (— l) n/ acquired by a vortex that encir- 
cles n f fermions. Because of this sensitivity, vortex tun- 
neling is potentially distinguishable on the background 
of more frequent fermion tunneling events. 

The contribution to tdc from fermion tunneling is sim- 
ply + 1)""£/,2> to lowest order in the fermion tunnel 
amplitudes tfi and tf 2 at the first and second constric- 
tion. There is no dependence on rif, so we need not 
consider it further. 

To calculate the contribution to tdc from vortex tun- 
neling, we apply the vortex tunnel Hamiltonian [24 
Hi = ViGi<j' i: where i = 1,2 labels the two constrictions 
and i>i is the tunnel coupling. The operators Oi and a[ 
create a vortex at the left and right end of constriction 
i, respectively. The lowest order contribution to td c is 
of second order in the tunnel Hamiltonian, because two 
vortices need to tunnel in order to transfer a single Ma- 
jorana fermion. The calculation of td c will be presented 
elsewhere, but the n v and rif dependence can be obtained 
without any calculation, as follows. 

Three terms can contribute to second order in Hi, de- 
pending on whether both vortices tunnel at constriction 
number 1 (amplitude tf), both at constriction number 2 
(amplitude t 2 ,), or one at constriction number 1 and the 
other at constriction number 2 (amplitude 2tit 2 )- The 
resulting expression for td c is 



t d c = t\ + t 2 + {-l) n f2t 1 t 2 , \in v 



is even. 



(13) 



We see that if the two constrictions are (nearly) identi- 
cal, so ti w t 2 = t, the tunnel amplitude tdc and hence 
the current /vortex due to vortex tunneling vanish if the 
fermion parity is odd, while / V ortox = (e 2 /h)V x At 2 if the 
fermion parity is even |26j . 

In summary, we have proposed a method to convert a 
charged Dirac fermion into a pair of neutral Majorana 
fermions, encoding the charge degree of freedom in the 
relative phase of the two Majorana's. The conversion can 
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be realized on the surface of a topological insulator at a 
junction between a magnetic domain wall (supporting 
a chiral charged mode) and two magnet-superconductor 
interfaces (each supporting a Majorana mode). We 
found that at low voltages the Dirac-Majorana conver- 
sion is geometry independent and fully determined by 
the electron-hole symmetry. It allows for the electrical 



read-out of a qubit encoded nonlocally in a pair of vor- 
tices, providing a building block for measurement based 
topological quantum computation. 
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the Dutch Science Foundation NWO/FOM. 
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